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We study the localization volumes V (participation ratio) of electronic wave functions in the 2d- 
Anderson model with diagonal disorder. Using a renormalization procedure, we show that at the 
band edges, i.e. for energies E ~ ±4, V is inversely proportional to the variance (e 2 ) of the site 
potentials. Using scaling arguments, we show that in the neighborhood of E — ±4, V scales as 



V = < £ 2 }- 1 !? ((4 
scaling ansatz. 



\E\)/(e }) with the scaling function g(x). Numerical simulations confirm this 



I. INTRODUCTION 

A large amount of work has been done in the past 
decades to understand the localization behavior in ran- 
domly disordered systems. The standard model for a 
single-particle electronic wave function in tight binding 
approximation in the presence of disorder is the Ander- 
son model [1-3]. In d = 2 and with diagonal disorder, 
it can be written as 

1pn+l 

= (E - 4)V»n,m - e n ,mVWn- (1) 

Here E is the energy, the hopping potentials between 
nearest-neighbor are all set to unity, (n, to) are the site 
indices, ifj is the eigenfunction, and |Vv?n| 2 is the prob- 
ability to find an electron at site (n, to). The e n . m are the 
site potentials which are all uncorrected random num- 
bers with the variance (e 2 ) = (1/iV 2 ) £^ m=1 e 2 n m . 

Their average value (e) = (1/iV 2 ) J2n,m=i e «,™ is set 
equal to zero. The term — 4^ n m on both sides of Eq. (1) 
is introduced in order to create a discretized Laplace op- 
erator on the left-hand side of the equation (see below). 

It has been recognized for long that in d = 1 and 
d = 2 all eigenstates of Eq. (1) are localized, whereas 
a localization-delocalization transition occurs in d = 3. 
However, the shape of the wavefunctions and the value 
of the localization length X(E, (e 2 )) is still being dis- 
cussed. 

In d = 1 and for uncorrelated site potentials, expo- 
nential localization was proven throughout the energy 
band [1,4,5] and a lot of rigorous results and scaling 
theories exist for the localization length A, defined via 
the Lyapunov exponent. Close to the band edges (i.e. at 
E = ±2 in d = 1), a weak disorder expansion yields 



X=(e 2 )- a f 



Ec - \E\ 



(2) 



with E c = 2, a = 1/3 and = 2/3 [6,7]. Recently, 
it has been shown by a space renormalization procedure 
[8,9] that Eq. (2) also holds for the case of long-range 
correlated site potentials with correlation exponent 7, 
< 7 < 1. In this case, the exponents have to be 



replaced by a — 1/(4 — 7) and (3 = 2/(4 — 7) an d 
7=1 refers to the uncorrelated case of Refs. [6,7], At 
the band center on the other hand, a different behavior 
of A occurs. A Green's function technique [10] yields 
A ~ (e 2 )" 1 for E = and in some distance from the 
band center, a second-order perturbation theory of the 
diagonal elements of the Green's function [1,11] yields 
\(E)~(4-E 2 )/(e 2 ). 

For d = 2, on the other hand, no analytical theory 
for the localization behavior is known yet. Numerical 
simulations close to the band center exist on the basis 
of Green's functions calculations [12], exact diagonal- 
ization [13] and the Lanczos algorithm [14], but do not 
lead to an exact or scaling form of A or related quanti- 
ties. Moreover, it was shown in Ref. [14], that the wave 
function in the two-dimensional Anderson model does 
not decay exponentially. Instead, a subexponential de- 
cay of ip was found with A increasing logarithmically 
with the distance r from the localization center. 

In this paper, we concentrate on the band edges, i.e. 
on energies E ps ±4 of the Anderson model in d = 2 
with uncorrelated potentials. We develop a renormal- 
ization approach, similar to the one in d = 1 of [8,9] 
and use it to find a scaling form for the localization vol- 
ume V, which is related to the inverse participation ratio 
P _1 . In d = 2 and with the wave function VVm be- 
ing normalized by J2 n ra ^ m = 1. P^ 1 is defined by 
[15,16] 



P- 



N 

E 



(3) 



Its inverse value P is a d-dimensional volume and mea- 
sures the extension of a given state. If we divide P by 
the volume Vq of the system, we get the relative vol- 
ume V of the eigenstate, V = P/Vq, i.e. the portion 
of the system where the wave amplitude is large. It can 
be easily verified that in d = 1,2 or 3, V ~ X d for 
all wavefunctions of the form ip(r) ~ cxp[— (r/X) ], 
$ > 0. Therefore, one can define an effective localiza- 
tion length V x l d ~ A, which measures the average di- 
ameter of the state. For numerical calculations in d = 2 
and d = 3, where the wave functions do not decay ex- 
ponentially, V is easier accessible than A and therefore, 



we focuse on V in this paper. 

Since V <~ A in d = 1, Eq. (2) holds up to an irrele- 
vant proportionality factor also for V. It is the purpose 
of this paper to show that a similar scaling law as the 
one of Eq. (2) holds also for V in d = 2, 

" -<«*>-•« (^j£) < 4 » 

but with different exponents, a = 1 and /? = 1 and with 
_E C = 4. This scaling ansatz is confirmed by numerical 
calculations. 

The paper is organized as follows: In section II we ex- 
plain the outline of the renormalization approach, while 
in section III the scaling ansatz for the localization vol- 
ume V is developed and tested by numerical simula- 
tions. Additional remarks about former renormalization 
theories and the extension to the vibrational problem are 
given in section IV. 

II. THE RENORMALIZATION APPROACH 

In the renormalization approach, we want to combine 
single sites of the lattice to blocks. This procedure must 
be reasonable in the limit of small values of 4 — \E\, 
i.e. close to the band edges. In this context, we must 
recall that the wavefunctions possess two characteristic 
lengthscales, (i) the wavelength A of the ordered lattice 
that describes the periodic fluctuating part of the wave- 
function and (ii) the localization length that describes 
their decaying envelope. As in d = 1 [9], we assume 
that the periodic part of ip does not depend on the disor- 
der but is reminiscent of the functions of an ordered lat- 
tice, where the disorder terms ejj are zero and ijj is a reg- 
ular sin- or cos-function. By solving Eq. (1) for the or- 
dered case, we found in d = 1 that A 2 <~ (\E C \ — E 1 ) -1 . 
In d = 2, where A 2 = A 2 + A 2 with the wavelengths A^ 
and A y in x— and y— direction, respectively, we have 

A 2 ~ (E c - IE])- 1 and A 2 ~ (E c - 

(5) 

where E c = 4. At the band edges, A diverges and the 
wavefunction no longer resolves the details of the disor- 
der potentials. In this case, we can imagine that neigh- 
boring sites of the lattice move as blocks and the follow- 
ing renormalization approach becomes legitimate. 

In the following, we consider the upper band edge 
E = 4, but by canceling the terms — 4^„ jTO on both 
sides of Eq. (1) and taking into account that the e n , m 
are randomly distributed around their mean value of 
(e) = 0, we can see that the equation is symmetric under 
the transformation E — > — E. Therefore, the renormal- 
ization approach is also valid for the lower band edge 
E = -4. 

In order to transform Eq. (1) into block form, we 
first replace the site indices (n, m) of the central site 
in Eq. (1) (see also Fig. 1) successively by (n + 1, m), 



(n — 1, m), (n, m + 1) and (n, m — 1). Combining those 
four equations with Eq. (1) and rearranging the terms, 
we arrive at 

(4/n,m ~t~ fn,m+l ~t~ fn,m—l H~ ,/n+l,m H~ fn—l,m) 
+ (E - 4)(4V>„,m + V'n+l.m + V>n-l,m + 

+4>n,m+l + VVm-l) + 8 V"n,m 
- 2(V>„+l, m +l + 1p n +l,m-l + V'n-l.m+l + V'n-l.m-l); 

(6) 

with the abbreviation /i j = ei,jipi,j. Comparing this 
result with Eq. (1), we can see that the left-hand side 
of Eq. (6) is again a Laplace operator, but with twice 
the distance between ip n . m and its neighbors. The first 
two terms on the right-hand side, involving the disorder 
terms foj and the eigenvalue (4 — E), are similar to the 
corresponding terms in Eq. (1), with the only difference 
that they no longer depend on a single site (n, m) but 
couple sites at distances < 2 from (n, m) to blocks. 

The last two terms, however, involve couplings be- 
tween ip n ,m and its second nearest neighbors ip n +i,m+i, 
V'n+i.m-i an d so on and do not occur in Eq. (1) (nor in 
the corresponding derivation in d = 1). Using a second- 
order Taylor expansion, we approximate these terms by 

^n+l,m— 1 ~l~ Ipn— l,m+l ~l~ VVi— l,m— 1 
~ -4V>n,m + 2(V> n+ i >m + VVi-l.m + 1pn,m+l + VVi.m-l) 
= 4?A„,m - 2/„, m + 2(E - 4)Vv™, (7) 

where Eq. (1) has been inserted in the last step. The 
Taylor expansion is legitimate in the limit of large wave- 
lengths, i.e. close to the band edge. Inserting Eq. (7) into 
Eq. (6) we finally arrive at 

(,/n,m+l H~ In , m — 1 H~ fn-\-l,m ~t~ fn— l,m) 

+ 1pn+l,m + 1pn-l,m)- (8) 

Assuming that the potentials are randomly distributed, 

(2) 

we introduce the smoothed wavefunction ip n ,m of the 
block and combine the terms f n , m +i + /n,m-i + 
/n+i,m + /n-i,m to one single term = e^m^m 

(2) 

with the block potential e n ' m = e n ,m+i + £n,m-i + 
en+i,m + £n-i,m- Equation (8) now shows a block form 
of block length ^ = 2, 

= -e^ m ^l + 4(E-4)^ m . (9) 

This is shown in Fig. 1 . Couplings between nearest- 
neighbor sites via Eq. (1) are symbolized by straight 
lines whereas the couplings between the sites ip n +2, m 
and ip n ,m and so on of Eq. (9) are symbolized by the 
oval lines. The site potentials / n , m+i = e„ iTO+i ^„ iTO+i 
and f n+ i tm = e n+iym ip n+iym with i = ±1 that form the 



block potential e„,,„ are indicated by the black circles. 
It can be seen that they lie well inside an inclined block, 
consisting of 2 2 particles. 
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FIG. 1.: Sketch of the Anderson lattice according to Eqs. (1) 
und (8) as explained in the text. The circles represent the dif- 
ferent lattice sites, the straight lines between them indicate 
the usual nearest-neighbor coupling, whereas the couplings of 
Eq. (8) are shown by the oval lines. The black circles stand for 
the site potentials that form the block potential [see Eq. (9)]. 
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FIG. 2.: Renormalization scheme for v — 4: the diagonal 
terms fij = ii,jipi,j that form the block potential are painted 
black and show a chess-board pattern. The sites (n,m), 
(n + 4,m), (n — 4, m), (n, m + 4) and (n, m — 4) that couple 
via a Laplace operator of distance v are symbolized by larger 
circles. 

This procedure can be continued. By replacing again 
the site indices (n, to) of Eq. (9) by (n + 1, to), (n — 
1, to), (n, ra+1) and (n, to — 1) and following the same 
procedure as before, we arrive at block indices v = 4. 
As long as the block length is well below A/2, the 
Taylor expansion is legitimate and we arrive at higher 
and higher orders of the renormalization. The poten- 
tial blocks form a chess-board pattern which is shown in 
Fig. 2 for the case ofiy — 4. The renormalized Anderson 
equation of block length v becomes 



(10) 



where V>£^m is the smoothed wavefunction of a block of 
length v and 



E 



t n,m — C *J 
i 2 +j 2 <v 2 A+j odd 



(ID 



with the sum running over all pairs of i and j with i even, 
j odd and vice versa (chess-board pattern) in a distance 
i 2 + j 2 < v 2 from the site index (n, to). 



III. THE SCALING ANSATZ: THEORY AND 
NUMERICAL SIMULATIONS 

Now, the renormalization approach is complete and 
we use it to derive a scaling theory for the localiza- 
tion volume V. Naturally, the form of the wavefunc- 
tion tp ~ exp[— (r/A)*] does not depend on the arbi- 
trary subdivision of the lattice into blocks. Nevertheless, 
by applying the renormalization approach over a certain 
range of block lengths, we gain information about V. 

The following derivation applies at E = 4, where A 
diverges and the block form is legitimate for any block 
size between 1 and infinity. At E = 4, the only quan- 
tities that enter into the right-hand side of Eq. (1) are 
the potentials e, ; j. Accordingly V, which is an average 
quantity over many lattice realizations, can only depend 
on the different moments of them (the first moment (e) 
being zero). As in d = 1 we presume a power-law be- 
havior, 



v~{e 2 y 



(12) 



To derive the exponent a, we apply the block trans- 
formation described above separately to both sides of 
Eq. (12). The left-hand side, V, is a volume and there- 
fore simply rescaled by a factor of v 2 , 



V -> V v 



V_ 

TT2 ■ 



(13) 



The right-hand side of Eq. (12) is determined by ran- 
dom walk theory. If we want to transform (e 2 ) into 
(e 2 )„, we must first summarize over all v 2 potentials 
Cij of one block and then calculate the variance over 
many different blocks. This is equivalent to calculating 
the mean square displacement of a random walk of v 2 
steps [17], 




v 2 (e 2 }. (14) 



Transforming Eq. (12) by Eqs. (13) and (14) we find 
V 



2a/,2\-a 



(15) 



As the last step, we must take into account that the block 
length v is arbitrary for A — > oo and Eq. (15) must there- 
fore not depend on v. This determines the exponent a 
and we finally find 



a 



1 and V-(e 2 )" 1 for E = 4. (16) 
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FIG. 3.: The effective localization length V 1/2 in the 2D- 
Anderson model at the band edge is plotted versus the vari- 
ance (e 2 ) of the site potentials in a double-logarithmic plot. V 
was calculated numerically for lattices of size 500 x 500 with 
Dirichlet boundary conditions and averaged over 40 systems. 
The line of slope —1/2 is a guide to the eye and shows the 
theoretical behavior. Finite size effects occur for small values 
of (e 2 ) (large values of V 1 ' 2 ). 



both quantities separately, but only on a suitable combi- 
nation of them. 

The scaling form of V can now be derived by standard 
techniques. Choosing v = (e 2 ) -1 / 2 , (which is permit- 
ted for large (e 2 ) even if A is not infinite) we find 



\E\ 



V{{±-\E\),{e 2 ))~{e 2 )^g 



with the scaling function g(x) and the argument 



x = 



4- \E\ 



(18) 



(19) 



For \E\ = 4, Eq. (18) must reduce to Eq. (16), yield- 
ing g(0) = 1 for x = 0. For small values of 4 — \E\ 
(large A) or large values of (e 2 ) (small V), VV <C A 
and x <C 1. In this case, the effective localization length 
\fV is smaller than A, the system behaves as if A were 
infinite and g(x) should therefore be a constant func- 
tion. For x 3> 1, the maximum block size becomes 
smaller and smaller, so that gradually, the scaling theory 
must break down. However, as in d = 1, an intermedi- 
ate range may exist, where g{x) still shows a power-law 
behavior. 



10 J 



In order to test Eq. (16), the eigenfunctions of sys- 
tems of size 500 x 500 with varying variances (e 2 ) have 
been calculated by the Lanczos algorithm. The differ- 
ent V have been determined using Eq. (3). For each 
(e 2 ), we took the average over 40 systems and calcu- 
lated the eigenfunctions in a small energy interval of 
E = 4 ± 0.0002. The results are shown in Fig. 3, where 
V is plotted versus (e 2 ) in a double-logarithmic way. 
The line of slope —1/2 is a guide to the eye and rep- 
resents the result of the scaling theory [see Eq. (16)]. 
Apart from slight finite-size effects for small (e 2 ) (and 
therefore large V 1 / 2 ) it agrees very well with the numer- 
ical results. 

The scaling theory can be extended to energies in 
some (small) distance from the band edge, where A 
is still large enough to perform the renormalization 
scheme over many steps. In Eq. (10), (4 — E) is rescaled 
with v 2 . Accordingly, we have (cf. Eqs. (13) and (14)) 



K=i {{±-\E\),{e 2 )) ~ v 2 V v {v 2 {l-\E\ 



(17) 



Equation (17) is a generalized homogeneity relation. 
This means that the form of V remains unchanged when 
both, (e 2 ) and 4 - \E\, are rescaled according to the 
renormalization theory. Thus, V does not depend on 
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FIG. 4.: As a test of Eq. (18), (V(e 2 )) 1/2 is plotted versus 
the argument x = (4 — E)/ (e 2 ) for different disorder widths 
w = 0.6 (circles), w = 0.8 (squares), w = 1.0 (diamonds), 
w = 1.2 (triangles up) and w = 1.5 (triangles down) with 
(e 2 ) = to 2 /12 and for different values for 4 — E between 
0.001 and 0.06. The average was again taken over 40 systems 
of size 500 x 500. 

In order to test Eq. (18), we have plotted (V (e 2 )) 1/2 
as a function of (4 — E)/(e 2 ) for different disorder 
widths w of the potentials, ej G [— w/2,w/2] with 
(e 2 ) = w 2 /l2 and for different values of 4- E. The nu- 
merical simulations were again carried out on 500 x 500 
lattices and the average was performed over 40 systems 
and in an energy interval [4 - £- 0.0002, 4 - E+ 0. 0002] 



for different values of 4 — E between 0.001 and 0.06. 
With decreasing values of (e 2 ), V increases and finite 
size effects occur. Additionally, systems where (e 2 ) and 
(4 — E) are both small or both large, possess large er- 
ror bars, i.e. large fluctuations between different values 
of V. In the case of small (e 2 ), this also gives rise to 
finite size effects, because some very large values of V 
are suppressed by the finite system size. So, we restrict 
ourselves to not too large values of 4 — E and to combi- 
nations, where such large fluctuations do not occur. 

The results are shown in Fig. 4 and confirm the scal- 
ing ansatz (18) very well. Different symbols that indi- 
cate different (e 2 ) fall onto the same universal curve. We 
can see that indeed reaches a plateau, g(x) w const 
for small values of x where V is simply described by 
Eq. (16) (see above). For large values x ^> 1, on the 
other hand, the scaling theory must break down, possi- 
bly after an intermediate range with a different power- 
law behavior of g(x). It would be very interesting to 
investigate also this regime, However, large values of x 
have not been calculated, because - due to the increas- 
ing vaues of \fV - we needed much larger system sizes 
for the simulations. This is currently not possible. 



with the eigenfrequency u of the vibration, we find the 
vibrational equation 

EOn'.m' - i>n,m) = -^VVm, (21) 




with the sum going over all neighbors of the site (n, m). 
Inserting the above transformation into Eq. (5), we 
find for the wavelength in the vibrational case A ~ 
((m) 1 / 2 ^) -1 . The limit of long wavelengths applies 
thus for uo 2 < l/(m). Positive masses lead to (to 2 ) < 
(m) 2 and together with Eq. (19) we finally arrive at 
x > 1. 

So, in the vibrational case, only the branch of higher 
values of the scaling variable x exists and it will be very 
interesting to investigate also this part. However, since 
this demands much larger system sizes (due to the in- 
creasing values of V), this should be done in the future. 
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IV. CONCLUSIVE REMARKS 

In summary, a renormalization scheme has been de- 
veloped close to the band edges (i.e. in the limit of 
large wavelengths) that analytically reduces the Ander- 
son equations (1) into block form where the block sizes 
may become arbitrarily large at the band edges. A scal- 
ing form for the localization volume V has been derived 
from this. Contrary to former renormalization schemes 
[18,19], it does not involve sucessive recalculations of 
the matrix elements in each step, but simply replaces 
Eq. (1) by Eq. (10), where the off-diagonal elements are 
unchanged and the diagonal elements of arbitrary block 
size v are directly related to the diagonal elements of 
the original system. The works of [18,19] proposed the 
mobility edge in d = 3 and scaling laws for the con- 
ductivity and related quantities. Therefore, it will be 
very interesting to extend also the present theory to the 
three-dimensional Anderson model. However, as it is 
developed for energies close to the band edge it is for 
the moment not clear, if it can be applied to the vicin- 
ity of the mobility edge, where comparisions to former 
renormalization theories can be made. 

As a last remark, we would like to note that the 
regime x > 1 is also relevant to the vibrational prob- 
lem with unit spring constants and fluctuating masses 
= (m) + fn n ^ m , where (m) describes the average 
mass and m niTO the disorder of them. If we transform 
Eq. (1) according to 

4 - E — > (m)uj 2 , e n , TO — ► fh n ^ m uj 2 (20) 
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